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TECHNICAL NOTE 4115 


THEORY OF AIRCRAFT STRUCTURAL MODEIS SUBJECT TO 
AERODYNAMIC HEATING AND EXTERNAL LOADS 1 
By William J. O'Sullivan, Jr. 


SUMMARY 


The problem of investigating the simultaneous effects of transient 
aerodynamic heating and external loads on aircraft structures for the 
purpose of deter mining the ability of the structure to withstand flight 
to supersonic speeds is studied. By dimensional analyses it is shown 
that: 


(a) A structural model geometrically similar to the aircraft and 
constructed of the same materials as the aircraft will be thermally sim- 
ilar to the aircraft with respect to the flow of heat through the structure 

(b) The thermal stresses and deformations of the structural model 
will be similar to those of the aircraft when the structural model is 
constructed at the same temperature as the aircraft. 

(c) The stresses and deformations of the structural model due to 
external loads will be similar to those of the aircraft. 

(d) By aerodynamic means the structural model can automatically be 
subjected to heating and cooling that correctly simulate the aerodynamic 
heating of the aircraft, except with respect to angular velocities and 
ang ular accelerations, without requiring determination of the heat flux 
at each point on the surface and its variation with time. 

(e) The similitude of aerodynamic forces, moments, and pressures 
acting on the aerodynamically heated structural model to those acting 
on the aircraft is determined for the case of zero angular velocity and 
zero angular acceleration, so that the structural model may be subjected 
to the external loads required for simultaneous simulation of stresses 
and deformations due to external loads . 
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HJTRODUDT ION 


The phenomenon of - aerodynamic heating of supersonic airplanes and 
guided missiles imposes the new design requirement that every component 
of the aircraft be capable of withstanding the adverse effects of the 
elevated temperatures and rates of change of temperature to which it is 
subjected in flight. Consideration is restricted herein to this problem 
with respect to the aircraft structure . 

Structural materials exhibit markedly different variations of 
strength with increase in temperature. Below the temperature of complete 
strength loss, the aircraft structure may experience serious decrease in 
strength and suffer other detrimental effects by a number of processes. 

If the structure were to have a uniform high temperature throughout but 
be composed of materials of different coefficients ■ of thermal expansion, 
stresses caused by unequal expansion would be generated in the structure} 
however, even if all parts had the same coefficient of thermal expansion, 
the rapid changes of velocity, altitude, and attitude of supersonic air- 
craft would produce transient aerodynamic heating which, together with 
finite thermal diffusivity (even of metallic structures), would result 
in unequal temperatures that would produce stresses as a result of unequal 
thermal expansion. The loss of material strength at elevated temperatures 
and the development of thermal stresses are attended by thermal distortion 
and altered elastic properties of the structure. These changes give rise 
to deflections that alter the externally applied aerodynamic loads as well 
as the distribution of load between the component parts of the structure. 
Dependent upon the type of structure, various kinds of thermally induced 
aeroelastic problems, such as panel flutter may result. 

If reliable computations of the transfer of heat between the air and 
every element of external surface of the aircraft at every instant of time 
during a flight and, in addition, computations of the flow of heat through 
the complex aircraft structure were possible, computations of the thermal 
stresses in the structure, the remaining strength, the thermal distortion, 
the alteration of aerodynamic characteristics and loads by the thermal 
distortion, the modification of the aeroelastic characteristics, and other 
important effects of transient aerodynamic heating might then be possible. 
Because of the enormousness of such a task, it appears that some other 
method of approach may be more practical. 

It is theoretically conceivable that by means of a dimensional anal- 
ysis the necessary scaling parameters might be developed whereby a struc- 
tural model of the actual aircraft could be built and subjected to 
1 aerodynamic- heating tests with attendant benefits analogous to those 
obtained with models in both free -flight and wind-tunnel aerodynamic 
investigations. Appropriate measurements upon such a structural model 
would yield the information sought, the physical phenomena involved being 
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thereby compelled in a sense to compute the wanted information. The pres- 
ent paper constitutes a first endeavor to develop rigorously by dimensional 
an a lysis the theory and attendant scaling laws for the construction and 
testing of such aerodynamic-heating structural models hav ing the requisite 
similarity to any given aircraft flying any prescribed path. 

For logical development of the theory, it has been found desirable 
to subdivide the presentation into the following parts: 

I. Thermal Similitude - This part includes the development of the 

equations defining the conditions under which there exists 
similitude of temperature, temperature distribution, heat con- 
tent, and flow of heat in the model structure to that in the 
aircraft structure. 

II. Thermal Stress and Deformation Similitude - This part presents 

the development of the equations expressing the conditions 
under which there exists similitude of thermally induced 
stresses and deformations between the model and aircraft struc- 
tures while thermal similitude is preserved. 

III. External Load Stress and Deformation Similitude - This part 
presents the development of the equations expressing the condi- 
tions under which there exists similitude of stresses and 
deformations due to external loads between the model and air- 
craft structures while thermal similitude and similitude of 
thermally induced stresses and deformations are preserved. 

IV. Heating the Structural Model - This part presents development 

of the equations expressing the conditions under which the 
structural model can be subjected, by aerodynamic means, to 
the requisite heating and cooling in simulation of the aero- 
dynamic heating experienced by the aircraft in flight, without 
requiring determination of the heat transfer to each point on 
the surface of the aircraft and model. 

V. External Loading of Structural Model When Aerodynamically Heated - 

This part presents the development of the equations relating 
the aerodynamic forces, moments, and pressures imposed on the 
structural model when aerodynamically heated to those of the 
aircraft to permit adjustment of the loads to loads required 
for similitude of stresses and deformations due to external 
loads. 
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SYMBOLS 

coefficient of thermal expansion of structure material, 
l/°R 

free -stream sound velocity, ft /sec 

exponents 

exponents 

specific heat of structure material, 

(Btu) (ft)/(lb) (sec2) (°R) 

free-stream specific heat at constant pressure, 
f t^/ ( sec2) ( °R) 

exponents 

mass density of structure material, lb-sec2/ft^ 
exponents 

modulus of elasticity in tension or compression, lb/ft 2 
exponents 

external concentrated load, lb 
modulus of elasticity in shear, lb/ft ^ 
heat, Btu 

aerodynamic heat-transfer coefficient, Btu/(ft^) (sec) (°R) 

aerodynamic moment, lb-ft 

thermal conductivity of structure material, 

Btu/(ft) (sec) (°R) 

free-stream thermal conductivity, lb/(sec) (°R) 

dimensional concept length 

linear dimension indicative of size, ft 
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dimensional concept mass 
external moment, ft- lb 
scale factor 

external distributed load, lb/ft 2 
aerodynamic pressure, lb/ft 2 
free-stream pressure, lb/ft 2 

dimensional concept heat 
heat flux, lb/ft-sec 
aerodynamic force, lb 
arbitrary exponent 
dimensional concept time 
time, sec 

free-stream velocity, ft/s ec 

angular orientation, radians 

angular velocity, radians/sec 

angular acceleration, radians/sec 2 

ratio of free-stream specific heats 

temperature change from assembly temper attire, °R 

thermal linear deformation, ft 

load linear deformation, ft 

dimensional concept temperature 

temperature, °R 

thermal diffusivity of structure material, ft 2 /sec 
free-stream absolute viscosity, lb-sec/ft 2 
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v Poisson's ratio 

free-stream mass density of gas, lb-sec 2 /ft^ 

or^ load tensile or compressive stress, lb/ft 2 

crp proportional limit in tension or compression, lb/ft^ 

(Tq thermal tensile or compressive stress, lb /ft 2 

t-q load shear stress, lb/ft^ 

Tp proportional limit in shear, lb/ft 2 

Tq thermal shear stress, lb/ft 2 

Subscripts: 

a aircraft 

m model 


I. THERMAL SIMILITUDE 


Formulation of Problem 

The objective of this section is to investigate what thermal sim- 
ilarity can exist between an aircraft structure and a model thereof under 
transient heating or cooling. By thermal similarity is meant the rela- 
tions of the temperature, temperature distribution, heat content, and 
flow of heat in the model to those in the aircraft structure. 


Physical Quantities Involved 

The assumption is made that the flow of heat inside the structure 
is predominantly by the process of thermal conduction through the solid 
material so that transfer of heat by radiation from one part to another 
and convective transfer of heat through the air inside the structure 
from one part to another are negligible. Accordingly, the transient flow 
of heat is assumed to be dependent upon the following physical quantities 
having the indicated consistent engineering _units and dimensions in terns 
of length L, mass M, time T, temperature 0, and heat Q. The con- 
sideration of temperature 0 and heat Q as distinct dimensions is 
justified in that in the phenomena here considered they retain their 
identity. 
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Symbol 

Quantity 

Units 

Dimensions 

l 

Linear dimension 
indicative of size 

ft 

L 

0 

Temperature 

°R 

0 

D 

Mass density 

lb-sec^/ft^ 

L _5 M 

C 

Specific heat 
of material 

(Btu) (ft) /(lb) (sec^) (°R) 


K 

Thermal conductivity 

Btu/ (ft) (sec) (°R) 


H 

Heat 

Btu 

Q. 

t 

Time 

sec 

T 


Dimensional Analysis 

By dimensional analysis (refs. 1 and 2), the conditions of transient 
thermal similitude are found as follows. Selecting length 1 , time t , 
temperature 0, mass density D, and heat H as the primary physical 
quantities, the resulting it-factors are 


0 I a l t b l e°l D d l B*1 
K l aa t' 32 9=2 D=2 3=2 

The exponents a^, a^, b-j_, b 2 > • • • are to be determined by the 

conditional equations upon L, M, T, 0, and Q. 

The conditional equations of the first it- factor containing C are 
obtained by substituting the dimensions of the physical quantities giving 

(m-%-^) (L) a l(T) bl (9) °l(l-3n) dl (Q) ei 


The conditional equations are then: 

L: a^_ - 5d^_ = 0 


M: 


-1 + dq = 0 
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T: b 1 = 0 

I 

0: -1 + Cy = 0 

Q: 1 + e-^ = 0 

from whence 

&! = 5 

t>i = 0 

c l= 1 
d!-= 1 

e ! = "I 

The dimensionless parameter resulting from the first rt-factor is 
then C130D/H. 

Similarly, substitution of the dimensions into the second jt-factor 
containing K gives 

(if W-^CL) *2(1) b 2(0) C2 ^L-5M) d2 (Q) e 2 

The conditional equations are then: 

L: -1 + ag - Jd 2 = 0 

M: d 2 = 0 

T: -1 + b2 = 0 

0: -1 + c 2 = 0 

Q: 1 + e 2 = 0 
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from -whence 


&2 = 1 
b2 = 1 
c 2 = 1 
d-2 = 0 
e 2 = -1 


so that -the dimensionless parameter resulting from the second it-factor 
is Klt0/H. 

The solution yielded by the dimensional analysis is then 



Klte\ 

E J 


0 


CD 


Similitude Equations 

For purpose of identification, the aircraft structural quantity is 
denoted with the subscript a and the model structural quantity with the 
subscript mj the similitude equations relating the transient heating of 
the aircraft and model structures may then immediately be written from 
equation (l) as 


Kja^m4m®m _ ^a^a^a^a 
= ^ 


( 2 ) 

(5) 


Equations (2) and (5) must be simultaneously satisfied. The condi- 
tions under which this requirement is possible may be examined as follows. 
Substituting equation (2) into equation (5) yields 


Kin'^m K a t a 
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Consider that, for every point in the aircraft structure, there must be 
a corresponding point in the model structure; then, a scale n between 
the aircraft and model structures may be defined as * 


Z m = n 2 a 


(5) 


Substituting equation ( 5 ) into equation (4) gives 


Km 

CmRn 


tm — 


n2 


J^-t a 

CaPa 


(6) 


In equation (6) it is observed that transient thermal similitude is 
dependent upon the ratio of thermal properties of the materials repre- 
sented by K/CD. This ratio is recognized’ as being the thermal diffu- 
sivity k so that equation (6) may be written as 


K m"km - n ^ K a^a 


( 7 ) 


In the analytical theory of the transient conduction of heat in solids 
(refs. 3 and 4), it is shown that the phenomenon is dependent only upon 
the one physical property of the material of thermal diffusivlty. Thus, 
this result indicated by equations (6) and ( 7 ) might have been foreseen. 
The thermal diffusivlty k , or its components, particularly the thermal 
conductivity K and the specific heat C, have received little investi- 
gation especially for structural materials. This is probably due in 
large measure to the fact that such investigations as have been made, 
primarily in the fields of calorimetry and metallurgy, have revealed that 
most materials exhibit significant anomalies in these thermal properties 
as functions of temperature. Consequently, it appears impractical to 
fulfi ll equation ( 7 ) over any appreciable range of temperature other than 
by constructing the model of identically the same material as used in the 
aircraft. Assume, then, that corresponding parts of the model and air- 
craft are constructed of identical materials, so that 


% = K a 

(8) 

Cm = Ca 

(9) 

4*m = E a 

( 10 ) 
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or 


K m - K a (ll) 

which implies, since these thermal properties are assumed to be functions 
of temperature, that 

®m = ®a (l2) 

The substitution of equations (8), ( 9 ) , and (10) into equation (6); or 
the substitution of equation (ll) into equation ( 7 ) then gives 

t m = n2t a (13) 

This means that an interval of time t a on the aircraft corresponds to 
the interval of time n^t a on the model. Thus, if the model is half 
the size of the aircraft so that n = anything that happens on the 

aircraft in one second will happen on the model in 0.25 second. The sub- 
stitution of equations (8), ( 9 ) , (lO), (12), and ( 13 ) into either equa- 
tion (2) or ( 3 ) gives’ 

Hm = n5H a (lit) 

which states that the heat content of the model iB v? that of the air- 
craft. This result is to be expected because the heat capacities of 
geometrically similar bodies composed of the same material and at the 
same temperature are proportional to their masses, which in turn are 
proportional to the cube of their sizes. 


Effects of Dissimilar Materials 

The foregoing analysis haB presumed the structure to be composed of 
the same material throughout. If different materials are used, then the 
physical properties of each material in the structure must be included 
R.mnng the physical quantities involved. One material can be selected 
as a reference material, and use will then be made of dimensionless param- 
eters that express the ratio of the physical properties of the reference 
material to the corresponding physical properties of each of the other 
materials . These requirements are automatically fulfilled when each part 
of the model is made of the same material as the corresponding part of 
the aircraft structure. The foregoing equations may therefore be con- 
sidered to apply to structures composed of any number of dissi m ilar 
materials . 
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Summary of Conditions of Thermal Similitude 

Under transient heating or cooling, thermal similitude •will exist 
between an aircraft structure and a model thereof when (eqs. (5), (8), 
(9) , and (10)) 

l m = nl a 
Km = K a 

Cm = C a 
Dm = D a 

The thermal similitude will then be (see eqs. (12), ( 15 ), (l4)) 

®m = ®a 

•tm = n2 t a 

= n ^Ha 

Equations (8), ( 9 ) , and (10) are automatically fulfilled when the model 
structure is made of the same material as the aircraft structure. 

From these similitude equations other useful similitude equations 
may be derived. Thus, from equations (12) and ( 5 ) , respectively, 

A8 m = AS a 


and 


Aim — n AZ a 


Dividing the first by the second and taking the limit yields 


d0_\ = i/de\ 

n W a 


(15) 


Thus, the thermal gradient in the model is l/n times that at the cor- 
responding point in the aircraft structure. . The similitude of the second 
derivative of temperature with respect to distance may be obtained as 
follows: 



Aim n A l a 
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In the limit. 



(16) 


Thus, the second derivative of temperature with respect to distance in 
the model is equal to l/n^ times the second derivative of temperature 
with respect to distance at the corresponding point in the aircraft 
structure . 


II. THERMAL STRESS AND DEFORMATION SIMILITUDE 


Formulation of Problem 

The objective of part II is to investigate the similarity of thermally 
induced stress and deformation that can exist between an aircraft structure 
and a model under transient heating or cooling when the conditions requi- 
site for thermal similitude developed in part I are fulfilled. 


Physical Quantities Involved 

At the temperature at which an aircraft structure is fabricated, 
all the parts fit together so that the structure is free of stress when 
it is not carrying externally applied loads, which include its weight. 

If the structure is made of a single isotropic material, and its tempera- 
ture throughout is changed a uniform amount from the assembly temperature, 
every element of material within the structure will undergo by thermal 
expansion an identical change in dimensions. The elements will fit 
together as they did before the temperature change, and therefore no 
thermally induced stress is developed. Again, if instead of a uniform 
temperature change, the structure is subjected to a constant temperature 
gradient throughout, and if the coefficient of thermal expansion is assumed 
to be constant, each of the elements comprising the structure will undergo 
a thermal expansion such that all elements fit together, and this result 
indicates that no thermally induced stress has been produced. On the 
other hand, if a nonuniform thermal gradient is imposed, the expanded 
elements do not fit together, so that thermally induced stress is indi- 
cated to have been produced. Thus, for the structure composed of the 
same material throughout, which is the case treated in the analytical 
theory of thermal stress (ref. 5), it is apparent that thermal stress 
can develop only when the temperature distribution through the structure 
is other than constant or linear. 
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The foregoing argument may be extended to the more general, and 
more important, case of the structure composed of any number of dissimilar 
but isotropic materials as follows. If such a structure is subjected to 
a uniform change of temperature from its assembly temperature, the dis- 
similar elements will undergo dissimilar thermal expansion so that they 
no longer fit together; consequently, the development of thermally induced 
stress is indicated. If the structure is subjected to a constant tempera- 
ture gradient, the unequal expansion of dissimilar elements will vary 
with position in the structure. The imposition of a nonuniform thermal 
gradient produces not only unequal thermal expansion between dissimilar 
elements but also, as in the case of the structure composed of the same 
material throughout, expansion of similar elements such that they will 
no longer fit together. 

These considerations show that, in the general case of the structure 
composed of dissimilar materials, the development of thermally induced 
stress is in general dependent upon the departure of the temperature from 
the assembly temperature dB, the temperature gradient de/d l, the sec- 
ond derivative of the temperature with respect to distance d 2 e/dZ 2 , and 
for each material involved its coefficient of thermal expansion A and 
its elastic properties. The elastic properties are the modulus of elas- 
ticity in tension or compression E (Young's modulus), the modulus of 
elasticity in shear G (modulus of rigidity), and Poisson's ratio v. 

Only two of these three elastic properties need be considered in that 
they are theoretically related (ref. 5 ) by 



2(1 + v) 


For simplicity, as well as for the purpose of indicating in subsequent 
sections that all three of these elastic properties must be identical 
in corresponding parts of the model and aircraft structures, all three 
of these elastic properties will be retained. 

It is of practical importance to know the closeness with which the 
stresses approach the elastic limits. Accordingly, the proportional 
limit in tension or compression bp and in shear Tp are included. 

Since it will be shown that the tensile, compressive, and shear stresses 
in the aircraft and model structures are equal when the similitude require- 
ments are fulfilled, and since these requirements demand that corresponding 
parts of the model and aircraft structures be composed of identical mate- 
rials, it may be alternatively considered that the closeness of approach 
to the elastic limits is implicit in other similitude equations. Again, 
in order to emphasize that the proportional limits are among the properties 
defining likeness of material, bp and Tp are explicitly included. 
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Lastly, among the physical quantities included are the length l, 
indicative of size and position, the thermally induced linear deforma- 
tion 6 q, the thermally induced tensile or compressive stress Cq, and 

the thermally induced shear stress tq. 

For brevity, the physical quantities pertaining to each different 
material in the structure are not written but rather each quantity is 
understood to pertain to each different material. Accordingly, the phys- 
ical quantities involved, their units in consistent engineering .measure. 


and their 
ture 0, 

dimensions in terms of length 
and heat Q are as follows: 

L, mass M, time 

T, temper a- 

Symbol 

Quantity 

Units 

Dimensions 

l 

Linear dimension indicative 
of size 

ft 

L 

5 e 

Thermal linear deformation 

ft 

L 

a e 

Thermal tensile or compressive 
stress 

lb /ft 2 

L-^MT" 2 

T0 

Thermal shear stress 

lb /ft 2 

L-3MJ-2 

a P 

Proportional limit in tension 
or compression 

lb /ft 2 

L-Imi;- 2 

T P 

Proportional limit in shear 

lb/ft 2 

L" 2 *®" 2 

E 

Modulus of elasticity in 
tension or compression 

lb /ft 2 

L" 2 *® -2 

G 

Modulus of elasticity in shear 

lb /ft 2 

L-lMT-2 

V 

Poisson's ratio 

Dimensionless 

Dimens ionle s s 

£S 

Temperature change from 
assembly temperature 

°R 

0 

d0/d l 

First derivative of 

temperature with respect to 
distance 

°R/ft 

L" 2 © 

d 2 e/di 2 

Second derivative of tempera- 
ture with respect to 
distance 

OR/ft 2 

L -2 © 

A 

Coefficient of thermal 
expans ion 

l/°R 

0" 1 
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Dimensional Analysis 

By dimensional analysis (refs. 1 and 2), the conditions for thermal- 
stress and thermal-deformation similitude under transient heating or 
cooling are found as follows. The above listed physical quantities 
involved in the phenomenon are observed not to involve the dimension 
heat Qj therefore, the physical quantities selected as primary quanti- 
ties need only contain the dimensions length. L, mass M, time T, and 
temperature ©. Accordingly, the quantities s_e_lected as primary quanti- 
ties are length 2, modulus of elasticity in tension or compression E, 
modulus of elasticity in shear G, and coefficient of thermal expansion A. 
The resulting rt-f actors in which the exponents of 2, E, G, A are 

denoted respectively by a, b, c, and d with numerical subscripts 
corresponding to the number of the rt-f actor are then as follows: 

8 e 2 ai E bl G Cl A dl 
<r e 2 a 2 E b 2 G 02 A d 2 

t 9 2 a 3 E b 3 g c 3 A d 3 

CTp 2 a ^ E b ^ G C ^ A d ^ 

T p 2^ E b 5 G C 5 A d 5 

v 2 a 6 E b 6 G c 6 A d 6 
£0 2®? E b T G c 7 a 3 ? 

2 a 8 E b 8 g c 8 A d 8 
d2 

2®9 E b 9 g c 9 

d2 2 - 

In order to illustrate the determination of the exponents ''by solu- 
tion of the conditional equations, the jt-factor containing Oq is 

employed as an example. This ^-factor is 

So b 0 Co do 
ff 9 1 2 E 2 G 2 a 2 
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Substituting for Oq , l, E, G, and A their dimensions in terms of L, 
M, T, and 0 gives 

(l - 1 *® -2 ) (L) a2 ^L"%T" 2 ) b2 (L“ ;i MI?" 2 ) C2 (0“ 1 ) d2 
The conditional equations are then: 

L: -1 + &2 - b2 - C2 = 0 

M: 1 + b2 + Cg = 0 

T: -2 - 2b 2 - 2 c2 = 0 

0: -d2 = 0 

from whence a 2 = 0, b 2 = -C 2 - 1, and dg = 0. The dimensionless param- 
eter resulting from, this jc-factor is then. 

g e Q- C2 

e 0 +l 
E * 

where C 2 is any arbitrarily chosen number including zero. Choosing 
C 2 as zero for simplicity gives 

fe 

E 

The exponents so determined for the nine ^-factors are 


ap = -1 

b ± = 0 

o 

ii 

r-l 

O 

d! = 0 

a2 = 0 

o' 

rv> 

ii 

i 

H 

Cp = 0 

d 2 = 0 

a^ = 0 

b 5 = -1 

Cj = 0 

o 

ii 

a^ = 0 

b^ = -1 

c^ = 0 

o 

ii 

ii 

o 

b 5 = -1 

Ctj = 0 

II 

o 

o 

ii 

of 

O' 

7\ 

II 

O 

o 

II 

VC 

o 

o 

II 

O 

II 

□r f 
—4 ( 

II 

O 

a 

— j ( 

II 

o 

&J = 1 

a 8 = i 

o 

II 

03 

n 

o 

II 

cO 

o 

rH 

II 

CO 

= 2 

b 9 = 0 

Cg = 0 

II 

H 
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The solution yielded by the dimensional analysis is then 



0 e T e 


E 


E 


X 

E 


E 


ASA, — IA, 
d l 



(17) 


Similitude Equations 


With the aircraft structural quantities denoted by the subscript a 
and the model structural quantities by the subscript m, the similitude 
equations are immediately written from equation ( 17 ) as 


^9,m _ ^8, a 
^m ^ a 

g 8,m _ cr 8,a 
Em E a 

T 9 ,m _ T 8 , a 
^m E a 

= °Pj a 

Em E a 

lEiS = ) 

Em ■ Ea I 


Vm - v a 
AGmAjj! = A0 a A a 

^nAn = (~r') ^a-^a 
m yl^/a 



(— \ lm 2 Am 

wL 



j 


( 18 ) 


The thermally induced stress and deformation equations ( 18 ) are to 
be satisfied, if possible, simultaneously with the thermal similitude 
equations ( 5 ), (8), ( 9 ), (10), (12), (15), (l4), ( 15 ), and (l6). That 
this is possible is shown by substituting into equations ( 18 ) the thermal 
similitude equation ( 5 ) defining scale 


^m = nZ a 
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and the following elastic-properties assumptions which are an extension, 
of the thermal similitude assumptions expressed by equations (8), (9), 
(10), and (ll) : 


Em — 

E a 

(19) 

v m = 

v a 

(20) 

°P,m = 

°p,a 

(21) 

T p,m = 

T p,a 

(22) 


-A-a 

(23) 


Equations (19) to (25) are to be understood as applying to corresponding 


elements in the model and aircraft, 
tions (l8) reduce to 

5 0,m = 
T 0,m = 

: 


With these substitutions equa- 


nS 0,a 

(24) 

°0,a 

(25) 

T 0,a 

( 26 ) 

£e a 

(27) 


and 



The latter two equations involving dG/dZ and d^9/dZ^ are identical 
with the thermal similitude equations (15) and (l6). It is to be observed 
that al 1 higher derivatives of temperature with respect to distance will 
similarly be compatible. Equation (27) is not only compatible with the 
thermal s imi litude equation (12) , but is actually a less general form of 
equation (12) . Thus, equations (24), (25) , and (26) are the only new 
equations resulting from the dimensional analysis of the similitude of 
ther mal stresses and deformations. Since they are in no way incompatible 
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with the thermal similitude equations, it is apparent that both thermal 
similitude and thermal stress and deformation similitude are simultaneously- 
attained in the same model structure when the model is geometrically sim- 
ilar to the aircraft structure and every component of the model is made 
of identically the same material as the corresponding part of the air- 
craft structure. 

Equations ( 25 ) and ( 26 ) show, respectively, that the thermally- 
induced tension (or compression) stress and shear stress in the model 
will be identical to those stresses in the aircraft structure at cor- 
responding points. Thus, by appropriate heating tests upon a model 
structure, the ability of an aircraft structure to withstand the adverse 
effects of aerodynamic heating or cooling can be determined. 

Equation (24) relates the linear thermal deformations of the model 
and the aircraft;. The linear thermal deformations of the model are n 
times those of the aircraft structure. Angular deformations in radian 
measure are expressed as the ratio of two linear deformations, and hence 
the similitude of angular deformations is implicit in equation (24) . It 
is convenient to have an explicit equation relating the thermally induced 
angular deformation of the model 00, m to that of the aircraft 0e ja . 

This relation is obtained from equation (24) by writing 

( 5 9,m) ;L n ( 5 0,a) 1 

^8, m 
00, a 

(S0,a) 2 ( & e,a) 2 

from whence 

08, m = 00, a ( 28 ) 

Thus, the thermally induced angular deformations of the model and air- 
craft structures are equal. 


( 5 0,m) 

l 2 

[ 5 0,a) 

Ka) 

'l ( 5 8 ,*)i 


Summary of Conditions of Thermal Stress 
and Deformation Similitude 

When the structural model fulfills the conditions of thermal simil- 
itude expressed by equations ( 5 ), ( 8 ), ( 9 ), ( 10 ), ( 12 ), (15), (l4), (15)> 
and (l6), it simultaneously fulfills the requirements for thermally 
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induced stress and deformation similitude when (see eqs. (19) to (25) 
and (27)) 

Em = 


a P,m " °P,a 
T p,m = T p,a 
An — Aq_ 
dBm. = dB a 


The s imi litude of thermally induced stress and deformation will then be 
(see eqs. (2I1-), (28), (25), and (26)) 

&8,m = nSe, a 

08, m = 08, a 
cr 8,m = cr 0,a 
T 0,m 3 T 0,a 

III. EXTERNAL LOAD STRESS AND DEFORMATION SIMILITUDE 


Formulation of Problem 

The objective of part HI is to investigate the similitude of 
stresses and deformations caused by externally applied loads that can 
exist between an aircraft structure and a model thereof when the model 
fulfills the requirements for thermal similitude and thermally induced 
stress similitude developed in parts I and II. 


Physical Quantities Involved 

An aircraft in curvilinear flight is subject not only to aerodynamic 
forces a.nH moments,’ including the thrust, but also to inertial forces and 
moments arising from the acceleration of gravity and linear and angular 
accelerations of the aircraft. For simplicity, the linear accelerations 
due to gravity and motion are customarily added vectorially so that the 
absolute linear accelerations are obtained. Since an aircraft is a 
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distended, body as distinct from a point mass, the absolute linear accelera- 
tions will, in general, vary from point to point within it. Also, its 
mass is irregularly distributed. Each element of the aircraft structure 
considered as a free body must be in equilibrium. In accordance with the 
principle of D'Alembert (ref. 6), the resultants of the externally applied 
forces and moments on the element are balanced by the inertial forces of 
the element due to its absolute linear accelerations and its angular 
accelerations, respectively. Thus, the portions of the externally applied 
forces and moments effective in producing stress within the element are 
the applied external forces and moments less the inertial forces and 
moments. Accordingly, in static tests of aircraft structures simulating 
flight, the absence of inertial forces and moments due to motion is cor- 
rected for by alteration of the externally applied forces representing 
aerodynamic forces and moments. 

If, then, the assumption is made that it is possible in a static 
test to apply to an aircraft structure external forces and moments so 
as to generate identically the same stresses and deformations as occur 
in flight, the problem resolves to a determination of the similitude of 
stresses and deformations that can exist between such a statically loaded 
aircraft structure and a model thereof when the model is constructed in 
accordance with the requirements for thermal similitude and thermal stress 
similitude developed in parts I and II. Thus viewed, the physical quanti- 
ties involved, their -units in consistent engineering measure, and their 
dimensions in terms of length L, mass M, and time T are as follows. 

The dimensions of temperature 0 and heat Q are emitted in that they 
are not involved. Also, the proportional limit in tension or compres- 
sion arp and the proportional limit in shear Tp are included among 

the quantities involved for the same reasons as in the thermal stress 
and deformation analysis of part II. 


Symbol 

Quantity 

Units 

Dimensions 

l 

Linear dimension indicative of size 

ft 

L 

&L 

Load linear deformation 

ft 

L 

F 

External concentrated load 

lb 

IME -2 

F ■ 

External distributed load 

lb /ft 2 

L“%E- 2 

N 

External moment — 

lb -ft 

L%T~ 2 


Load tensile or compressive stress 

lb /ft 2 

L - ^ -2 

t L 

■Load shear stress 

lb /ft 2 

L“^MC“^ 

a P • 

Proportional limit in tension or 
compression 

lb /ft 2 

L-^C" 2 
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Symbol 

Quantity 


Units 

Dimensions 

T P 

Proportional limit in 

shear 

lb/ft 2 

L _1 MT- 2 

E 

Modulus of elasticity 
sion or compression 

in ten- 

lb /ft 2 

l - 4mt -2 

G 

Modulus of elasticity 

in shear 

lb /ft 2 

L-%0- 2 

V 

Poisson's ratio 


Dimensionless 

D imens ionles s 


Dimensional Analysis 

By dimensional analyses (refs. 1 and 2), the conditions of simili- 
tude of the external loads and moments and the associated structural 
stresses and deformations are found as follows. Since the physical 
quantities involve only the dimensions length L, mass M, and time T, 
the physical quantities selected as primary quantities are length l, 
modulus of elasticity in tension or compression E, and modulus of elas- 
ticity in shear G. The resulting jt-factors in which the exponents of l , 
E, and G are denoted, respectively, by a, b, and c with numerical 
subscripts corresponding to the number of the jt-f actor are then as 
follows : 


8l l &1 E b l G C 1 


F E b 2 G°2 

ai bz cz 

Pl^E^G 5 


N A E b 4 G C ^ 


1*5 

b= 
E 5 

G c 5 

l a 6 

E b 6 

G c 6 

i *7 

E b T 

G c 7 

l a 8 

E b 8 

G c 8 

“9 e”? a 9 
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Determination of the values of the exponents in the 
in part II yields: 


a l = 

-1 

bl = 

0 

ci = 

a 2 ~ 

-2 

b 2 = 

-l 

c 2 = 

a 3 = 

0 

cr 

va 

n 

-l 

c 3 = 

a^ = 

-3 

b 4 = 

-l 

c 4 = 

a 5 = 

0 

b 5 = 

-l 

c 5 = 

a 6 = 

0 

b 6 " 

-1 

c 6 = 

a^ = 

0 

b ? = 

-l 

c ? = 

f» 

CD 

ii 

0 

b 8 = 

-l 

Cg = 

a 9 = 

0 

b 9 = 

0 

C 9 = 


same manner as 

0 

0 

0 

0 

0 

0 

0 

0 

0 


The solution yielded by the dimensional analysis is then 



t L «p 
E * E * 




(29) 


Similitude Equations 


Denoting, as in preceding sections, the aircraft structure by the 
subscript a and the model structure by the subscript m, the similitude 
equations are immediately ■written from equation (29) as 


= S L,a 


■m 


^m^m 


F a 

*a 2 E 



fin _ P a 
Fm E a 


(50) 


(continued on next page) 
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& 



N a 

Zm^Em 

2a^Ea 

°L,m 

a L,a 

Em 

E a 

T L,m 

T L ? a 

Em 

E a 

ffp,m _ 

°P,a 

Em 

E a 

Tp,m _ 

T p,a 

Em 

Ea 


(30) 

Cone. 


'm = v a 


The external load., moment, stress, and deformation similitude equa- 
tions ( 30 ) are to be satisfied, if possible, simultaneously vith the 
thermal similitude equations ( 5 ), ( 8 ), ( 9 )> ( 10 ), (l 2 ), (lj), (l4), ( 15 ), 
and (l6) and with the thermal stress and deformation similitude equa- 
tions ( 19 ) to (28) . That this is possible and the conditions under which 
it is possible are found by substituting into equations ( 30 ) the thermal 
similitude equations and the thermal stress and deformation similitude 
equations, whereupon equations ( 30 ) reduce to 


5l,m 

= nS L , a 

(31) 

Fm = 

**a 

(32) 

Pm ; 

= Pa 

(33) 

Nm = 

n^a 

(3t) 

°L,m 

= 0 L,a 

(35) 


* 
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T L,m “ T L,a 


(36) 


Thus, when the model is geometrically similar to the aircraft structure 
and constructed of the same materials, the externally applied forces, 
pressures, and moments (corrected for inertial forces and moments) that 
must be applied to the aircraft structure in a static test properly to 
simulate those obtained in flight are related to the externally applied 
forces, pressures, and moments that must be applied to the model struc- 
ture by equations (52), (55), and (54-), respectively. When this rela- 
tionship is obta in ed, the elastic linear deformations of model and air- 
craft are related by equation (51), the tensile or compressive stresses 
by equation (55), and the shear stresses by equation (5 6) • 


The relation of model and aircraft angular deformations under these 
externally applied loads may be derived from equation (51) in a manner 
similar to the derivation of equation (28). Denoting the load induced 
ang ular deformation of the model by 0L,m and that of the aircraft struc- 
ture by 0L,a> ( 

( 8 L,m) 1 n ( & L,a) 1 

0L,m _ 

0L,a ( 5 L,a)i 

( 8 L,a ) 2 


n l 5 L,a ) 2 
( 5 L, a )i 

( S L,a ) 2 


and 


0L,m = 0L, a (57) 

It is to be noted that the pressure, force, and moment similitude 
equations (55), (52), and (54) are those associated with geometrically 
similar bodies by virtue of the fact that the model and aircraft struc- 
tures f ulf ill the requirement 


lm = nla 

Inasmuch as force is the product of pressure acting over an area. 


*m _ p m z m 2 = n2p a l a 2 
p a p a^a^ p a^a^ 
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whence 

F m = n 2 F a 

Since moment is the product of force and distance, 

jjm _ Pm7m _ n ^a^a _ ^ 

N a Pa7a Pa^a 

whence 

N m = n^Na 


Summary of External Load and Stress Similitude 

When the structural model fulfills the conditions of thermal simili- 
tude expressed by equations (5)> (8), (9), (10), (12), ( 13 ) , (1*0 j ( 15 )> 
and (l6) , and the thermal stress similitude expressed by equations (19) 
to ( 28 ), the model can then simultaneously fulfill similitude with respect 
to the aircraft structure statically loaded to correctly simulate flight 
inertial and aerodynamic loading under the conditions (see eqs. ( 30 ) 
to (37)) 


^L,m “ n 8L,a 


Fm = n^F a 
?m = Pa 

N m = n5N a 

a L,m = a L,a 
T L,m = T L,a 
0L,m = 0L,a 
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IV. HEATING THE STRUCTURAL MODEL 


Formulation of Problem 

The problem considered in part IV is that of applying to each element 
of surface of the structural model the requisite heating and cooling to 
simulate correctly the aerodynamic heating -which the aircraft experiences 
in flight. The model is assumed to be constructed in accordance with the 
equations of part I so that similitude with respect to the flow of heat 
through the structure exists between the model and the aircraft, in 
accordance with the equations of part II so that there exists similitude 
between the model and aircraft with respect to the thermally induced 
stresses and deformations, and in accordance with the equations of part III 
so that there also exists similitude between the model and aircraft with 
respect to the stresses and deformations produced by external loads. The 
thermal character of the model has been fixed by these considerations so 
that the problem of heating and cooling the model must be compatible with 
these prior requirements. 

These prioh requirements with respect to the transfer of heat to 
or from the model are embodied in two equations. If the surface heat 
flux q is defined as the quantity of heat per unit of area per unit of 
time, the relation between the model and aircraft heat fluxes is obtained 
from equations ( 5 ), (l 3 )> and (l4) by writing 

- ^ _ n ^a _ 1 

^ ^m^fcm n^a^^ta n 


or 



( 58 ) 


Thus, the heat flux into or out of any elemental area of the model sur- 
face must be l/n times that of the corresponding element of surface of 
the aircraft at corresponding instants of time given by equation ( 13 ), 
which is 


t m = n 2 t a 

In the present state of development of theoretical and experimental 
boundary-layer information concerning aerodynamic heating, it is not pos^ 
Bible to compute reliably the heat flux to each point on the surface of 
complex bodies like airplanes or guided missiles for a specified trajec- 
tory. If the heat flux is determined experimentally for the particular 
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body traversing a particular trajectory, there remains the problem of 
devising apparatus capable of applying to the model in correct time 
sequence the proper heat flux which, in general, may vary from positive 
to negative, that is, from heating to cooling. It is conceivable that 
if the parameters governing the heat flux are properly scaled, the req- 
uisite heat flux applied to the model by aerodynamic heat ing may pos- 
sibly be determined by means of a suitable wind tunnel or by free flight 
of the model. If such is possible, both the problem of determining the 
heat flux and that of applying the heat flux are simultaneously solved; 
accordingly, this possibility is investigated in part IV. 


Physical Quantities Involved 

In the light of present experimental and theoretical knowledge, the 
physical quantities involved in the phenomenon of aerodynamic heating of 
aircraft in the region of immediate interest, where radiative heat exchange 
is negligible compared to forced convective heat exchange, are the fol- 
lowing quantities and have the indicated consistent engineering units 
and dimensions in terms of length L, mass M, time T, and temperature 0. 
The dimension heat Q is not here used because the phenomenon considered 


involves the conversion of energy between 
heat state. 

the mechanical 

state and the 

Symbol 

Quantity 

Units 

Dimensions 

q. 

Local heat flux 

lb/ft-sec 

MT-3 

8w 

Local surface temperature of 
body 

°R 

0 

000 

Free-stream temperature 

OR 

0 

Poo 

Free-stream pressure 

lb /ft 2 

L-^" 2 

Poo 

Free-stream mass density 

lb-se^/ft^ 

l~3m 

Vco 

Free-stream velocity 

ft/sec 

1 

H 

Poo 

Free-stream absolute viscosity 

lb-sec/ft 2 

L-lMT-l 

koo 

Free-stream thermal 
conductivity 

lb/ (sec) (°R) 

IMC-30-1 

c p 

Free-stream specific heat at 
constant pressure 

ft 2 / (sec 2 ) (°R) 

L 2 !" 2 ©" 1 

7 

Ratio of free-stream specific 
heats 

Dimens ionle s s 

Dimensionless 
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Symbol 

l 

a 

a 

a 


Quantity 

Body linear dimension 

Orientation of body relative 
to stream 

Angular velocity 

Angular acceleration 


Units 

ft 

radians 

radians/sec 

radians/sec^ 


Dimensions 

L 

Dimensionless 

T-l 

T-2 


Dimensional Analysis 

A dimensional analysis may now be performed, the primary physical 
quantities selected being the velocity the gas density , the ref- 

erence linear dimension l, and the thermal conductivity k^. The fol- 
lowing it -factors result in which the exponents of p ro , l, and 

are denoted respectively by a, b, c, and d with numerical subscripts 
corresponding to the number of the it-f actor: 


a-i bn Cn dn 

v a 2 b 2 ? c 2 v d 2 
-P M ^ ^00 


bz Cz Cbr 

0«> v ffl 5 P m 5 1 5 koo 5 


a 4 *4 ^4 t d 4 

v =o P M l K> 


a= be Ce de 

Poo Voo l Ko 


= v “ 6 « ‘ 6 l' 6 L, 46 

C P Voo p_ fr Koo 


av \> r 7 C7 d>7 

rv.'p, 'i ' ' 


a v /8 p „ b 8 l =8 ^8 

i V„“9 p b 9 I°9 ^ 


a vA 10 p„ b l° Z° 10 k> 0 
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Solution of the conditional equations yields the values of the exponents 
as : 


aq = -3 

H 

I 

II 

H 

,Q 

Ci = 0 

P 

11 

0 

&2 = -3 

b 2 = -1 

c 2 = -1 

d 2 = 1 

a 5 = -5 

b 3 = - 1 

c 3 - - 1 

d 3 - 1 

a^ = -2 

= -1 

C4 = 0 

d 4 = 0 

*5 = "I 

b 5 = 'I 

c 5 = “ 1 

d 5 = 0 

a 6 = ! 

b 6= 1 

c 6= 1 

d 6 = -1 

o 

n 

r 

by = 0 

C rj = 0 

dy = 0 

00 

II 

o 

bg= 0 

Cg = 0 

dg= 0 

yg 

ii 

* i 

H 

b 9 = 0 

c 9 = 1 

dy = 0 

a 10 = “ 2 

biQ = 0 

c 10 = 2 

d 10 = 0 


The solution yielded by the dimensional analysis is then: 




(39) 


In this derivation of equation (39 ) , the dimensionless parameters 
have been directly derived in terms of the free-stream parameters, whereas 
in the analytic boundary-layer theory the parameters governing aerodynamic 
heating are usually expressed in terms of the local conditions at the 
point under consideration. It is to be observed that equation (59) may 
be alternatively derived by dimensional analysis by first finding the 
dimensionless parameters in terms of the local conditions and then, by 
a second dimensional analysis, finding their transformation to free- 
stream conditions . 
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Identity of Dimensionless Parameters 

The dimensionless parameters occurring in equation (59) are not 
immediately recognizable as corresponding to the dimensionless parameters 
customarily employed in heat transfer and in aerodynamics. This is 
because the customary heat-transfer parameters were originally developed 
from consideration of an incompressible fluid, and the aerodynamic param- 
eters were separately developed without regard to heat transfer. In equa- 
tion (59) both heat transfer and aerodynamics are combined. The first 
parameter q/Voo^p^ corresponds to the Nusselt number hl/k of heat 

transfer, where h is the forced convective heat-transfer coefficient, 
for had h been employed instead of q only the first parameter of 
equation (59) would be altered, being changed to the Nusselt number. In 
the second parameter QyKo/^oo^Poo^ ^ ' t ^ ie third parameter 

the quantity VoPp^/Ko has the dimensions of temperature, and thus the 

second and third parameters together correspond to the thermal ratio 8 w /0 
of the analytic theory of aerodynamic heating (ref. 7) . The fourth param- 
eter Poo/Voo^p^, corresponds to the Mach number Voo/ao, where aoo is the 

velocity of sound in the free stream. Substituting into the fourth param- 

p aoo 2 

eter the equation for the velocity of sound in a gas, — = , gives 

Poo 7 

So^/tVoo 2 . The ratio of the specific heats 7, being a constant, may be 
dropped from the dimensionless parameter. The remainder is the square 
of the reciprocal of the Mach number. Since the reciprocal of a dimension 
less parameter, or any power or root of the parameter or its reciprocal 
may be used with equal validity, the fourth parameter corresponds to the 
Mach number . The fifth parameter P00/V00P l is the reciprocal of the 

Reynolds number Voop^ / Moo • The sixth parameter CpVoop^z/k^ corresponds 
to the Prandtl number Cpp^/k^, for has the dimensions of 

as may be seen from the fifth parameter. The seventh parameter, the ratio 
of the specific heats 7, is a property of the gas employed that is neces- 
sitated by the compressibility of the gas and the conversion of energy in 
the flow, between the mechanical and heat states. For brevity, the eighth 
parameter a is to be understood as representing the several angles 
defining the orientation relative to the stream, such as the angle of 
attack and the angle of yaw. In keeping with this abbreviation, the last 
two parameters, &Z /V^ and al 2 /V„2, are to be understood as representing 
the several flow angles imposed upon the distended body by its several 
angular velocities and accelerations. These parameters are thus recog- 
nized as being associated with the several stability derivatives. 
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Similitude Equations 


With the aircraft quantities denoted by the subscript a, and the 
model quantities by the subscript m, the similitude equations may be 
written from equation (59) 85 



(40) 
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Inasmuch as the model is constructed so as to satisfy the thermal 
si mi litude equations of part I, the thermal stress and deformation simili- 
tude equations of part II, and the equations of similitude of stress and 
deformation due to external loads of part III, there may be substituted 
into equations (40) equation (5)j which is 


^m - n ^a 


and equation ( 58 ) > which is 

Qm = 5 <3a 

Also, if the model is assumed to be constructed at the same temperature, 
at which the aircraft is constructed, the equation ( 27 ) 


^m= ^a 


may be rewritten as 


9m ~ ®a (41) 

Since equation (4l) applies to any point in the structure, it follows 
that 

9w,m = e w, a (^) 

Substituting equations ( 5 ), ( 58 ), and (42) into equations (40), and then 
solving equations (40) simultaneously insofar as possible results in 


4oo,m - a 
®°° ; m = 9co ; s 



(m» c p) ra ~ (l i oo c p) a 
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7m = 7 a 


£% - 



If, for the moment, the last three of these equations are disregarded, 
it is apparent by inspection that the remaining five equations can be 
compatible with the conditions imposed upon the model in parts I, II, 
and III only if the model is tested in the same gas in which the aircraft 
flies, namely air, for it is not possible to specify at will the properties 
of gases. On the assumption, then, that the model is tested in air, 


c p,m ~ c p, a 

and by the general gas equation 


and 



% e oo,m 



R g e «, 


where E g is the gas constant for air. With these relations, together 
with the first of equations (!«)) and equation ( 58 ), the five equations 
reduce to 

kco,m = ko=,a ^3) 

9co,m = 9oo, a (^0 
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, a P 

00 , m n ^ 00 , a 

(45) 

a>,m n *«>,a 

(46) 

V = V 

v °°,m v oo,a 

(47) 

u = u 

^00 ^ a 

(48) 

c p,m = c p, a 

(49) 

11 

P 1 2 

(50) 


Consider now the three equations that were momentarily disregarded, 
which, inasmuch as Voo^ m = Voo,a.> ncrvr reduce to 

%= % (51) 

= ^ ^a (52) 

Sm - % (53) 

vr 

The first of these equations, = a^, is in no way incompatible, for 

it is physically possible at corresponding instants of time, as given 
by equation (13) , 


tm = n2 t a 

for the model and the aircraft to have the same attitude; but if (% = {% 
and % = n 2 t a , it follows by differentiation that 

1 • 

°m - ~o a u 
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Thus, the aerodynamic-heating similitude equation involving a is com- 
patible with the similitude equations of parts I, II, and HI, but the 
aerodynamic-heating similitude equations involving d and a are not 
compatible except when the scale factor n is unity. 

The significance of this incompatibility with respect to a and a, 
is that, for a model for which n ^ 1, it is not possible to simulate by 
aerodynamic means the effects upon the aerodynamic heating that are attrib- 
utable to angular velocities and angular accelerations. These effects are 
alterations of the direction of flow. For example, in regard to angular 
velocity, a fighter airplane, capable of a maximum rolling velocity of 
200° per second and having a wing span of 50 feet, when flying in the 
stratosphere could develop at its wing tip a change in angle of attack 
of about 2.6° at Mach number 2 and about I. 7 0 at Mach number 5- The 
other angular velocities and associated local flow-direction changes that 
the airplane could reasonably encounter at speeds at which aerodynamic 
heating is significant are of much lower magnitude. At the present state 
of development of theoretical and experimental information about aero- 
dynamic heating, it is not possible to evaluate the importance of such 
small flow-direction changes, but by virtue of their smallness their 
effects are probably small. The flow changes associated with angular 
accelerations are generally of an even smaller magnitude. 

Since the aerodynamic method of applying heating and cooling to the 
structural model is capable of exactly simulating the aerodynamic heating 
experienced by the aircraft in every respect except those of angular 
velocities and angular accelerations when n 1, which defect is of 
questionable significance even in the most violent maneuvers, it is 
apparent that -the method is highly desirable. 


Summary of Aerodynamic Heating S imi litude Equations 

A structural model constructed at the same temperature as the air- 
craft in accordance with the internal heat-flow similitude equations of 
part I, the thermally induced stress and deformation similitude equations 
of part II, and the external-load stress and deformation similitude equa- 
tions of part III can be subjected to aerodynamic heating and cooling 
exactly simulating those obtained in flight of the aircraft in every 
respect (except in regard to angular velocities and angular accelerations) 
when subjected to an airstream under the following conditions (see eqs. ( 13 ) , 
(44) to (47), and (51)): 

^m = n ^a 
®»,m = 9<»,a 
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■^oo, m ~ n -^ooj a 


P = — p 
K co,m n M ca,a 

V = V 
v oo^m v co ; a 

«m = “a 

The following conditions are then automatically satisfied (see eq_s. ( 38 ), 
(42), (43)., and (48) to ( 50 )): 

„ _ 1 „ 

‘Sm - n %. 

®w,m = ®w, a 

^00, m = ^00, a 

^,111 = a 

c p,m = c p,a 

?m = ?a 


Physical Possibility of Aerodynamic Heating Method 

The physical possibility of aerodynamically heating the structural 
model as distinct from the practical problems involved, which are beyond 
the scope of this investigation, are now considered. 

If the required time-varying air flow is to be supplied by a super- 
sonic wind tunnel, it is foreseeable that the density, pressure, temper- 
ature, and velocity in the test section may theoretically be obtained by 
extensions of the methods currently employed. Such a tunnel will nec- 
essarily provide a throat section for which the geometry can be varied 
during a run. This variable geometry might be accomplished by use of 
the sliding-block or flexible-wall type of adjustable supersonic wind- 
tunnel nozzle. 
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The possibility of employing a structural model in free flight is 
next considered. At any given time the atmosphere has a variation of 
density, pressure, and temperature with altitude that is not alterable 
at will. Equation ( h-h-) requires that the free-stream air temperature of 
the model and aircraft be the same, whereas equations (1*5) and (1*6) 
require that both the free-stream pressure and density of the model of 
scale n be l/n of that of the aircraft. These requirements can 
obviously be met for a model having a value of. n other than unity only 
in an isothermal atmosphere which is approximated in the stratosphere. 

If the standard aeronautical atmosphere (ref. 8) -and a half -scale model 
for which n = 0-5 are assumed, flight of the model at the tropopause 
(altitude, 56 , O 89 geopotential feet) corresponds to flight of the air- 
craft at an altitude of ^0 , 511 geopotential feet. It Is to be recognized, 
however, that many guided missiles are sufficiently small so that they 
may be tested as free- flight models at full scale and therefore would 
not require alteration of altitude from that at which the missile operates. 


V. EXTERNAL LOADING OF TEE STRUCTURAL MODEL 
WHEN AERODYNAMICALLY HEATED 


Formulation of the Problem 

The aerodynamic method of heating the structural model developed in 
part IV simultaneously applies aerodynamic loads to the model. A deter- 
mination of the similitude of these external aerodynamic loads imposed 
on the model to the aerodynamic loads encountered by the aircraft in 
flight Is therefore necessary in order that properly adjusted external 
loads may be imposed on the inode 1 for correct similitude of the stresses 
and deformations caused by external loads as discussed in part III. 


Physical Quantities Involved 


Frcm consideration of both theoretical and experimental information, 
the physical quantities influencing aerodynamic forces and having the 
in dicated consistent engineering units and dimensions in terms of 
length L, mass M, and time T are as follows: 


Symbol Quantity 

R Aerodynamic force 

l Linear dimension 


Units 

lb 

ft 

lb-sec^/ft^ 


Dimensions 

IMP - 2 

L 

L -5 M 


Free-stream mass density 
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Symbol 

Quantity 

Units 

Dimensions 

Voo 

Free-stream velocity 

ft/sec 

LT-1 

Poo 

Free-stream pressure 

lb/ft 2 

L -1 Mr' 2 

Moo 

Free-stream absolute viscosity 

lb-sec/ft 2 


koo 

Free-stream thermal 
conductivity 

lb/( sec) (°R) 

imt-3©- 1 

c p 

Free-stream specific heat at 
constant pressure 

ft 2 /(sec 2 ) (°R) 

L 2 ! -2 ©- 1 

7 

Ratio of free-stream specific 
heats 

Dimensionless 

D imens ionle s b 

9v • 

Local surface temperature of 
body 

°R 

0 

Qoo 

Free-stream temperature 

°R 

0 

a 

Orientation of body relative 
to stream 

. radians 

Dimensionless 

• 

a 

Angular velocity relative to 
stream 

radians/sec 

T-l 

• • 

a 

Angular acceleration relative 
to stream 

radians/sec 2 

T~2 


The inclusion of koo, Cp, 7 , 0 W , and 0oo is justified on the 

grounds that the viscous forces are influenced by the aerodynamic heating. 
The exclusion of the heat flux q. is in accord with equation (59), which 
shows that q is a function of the other quantities considered. For 
brevity, as before, a is to be understood as representing the several 
angles required to define the attitude of the body, and accordingly d 
and d are to be -understood as representing the first and second deriv- 
atives of these angles with respect to time. 


Dimensional Analysis 

A dimensional analysis may be performed by selecting as the primary 
physical quantities Vo,, p OT , l, and k*,. The following rt-factors 

result: 
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iT 


R v„ a i p w b i 2 C 1 k/l 

P. V » a 2 P> 1 C2 l ^ 2 
M» v/5 P J>3 i'3 ^3 
op V/4 p \ fb *> 
.°5 o b 5 l b 5 v ‘ *5 


7 V„ - p 


8 W v„ a 6 p b 6 l°6 kjt 


9. p„ b 7 l c 7 ka, d 7 
a V „ a 8 p J >8 j =8 k *6 


i C 9 p *” 9 1 ° 9 *.*» 


a O® P b i° i c i° t/io 


Solution of the conditional equations yields the values of the exponents 
in the % -factors as 


4? 

II 

to 

H 

1 

II 

H 

ja 

0 
H 

II 

1 

ro 

d. 1 = 0 

ag = -2 

b£ — -1 

o . 
n 

CM 

O 

o 

II 

-d 1 

= —1 

b 5 = -1 

c 3= ‘I 

d 5 = 0 

H 

II 

of 

\ - 1 

Cjj. = 1 

H 

1 

II 

a^ — 0 

b 5 = 0 

c^ = 0 

vf 

II 

o 

4f\ 

1 

II 

b 6- -1 

c 6 = -1 

H 

II 

-tP 

8sj — -5 

b^ = -1 

0 
-M 

II 

1 

H 

&y = 1 

p 

CD 

II 

O 

o 

ii 

CO 

O 

II 

CO 

o 

P 

CD 

II 

o 

v£ 

II 

1 

H 

b 9 = 0 

c^ = 1 

o 

II 

a 10 " “ 2 

o 

II 

.o 

c 10 = 2 

d 10 - 0 
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The solution yielded by the dimensional analysis is then 





V_ 2 p 

00 * c< 




K> 


7 , 




vjpl vjpl 



~r ] = 0 (54) 


It is observed that, with the exception of the first dimensionless 
parameter of equation (54) , the parameters are the same as those of equa- 
tion (59) (but for the absence of the q/Voo^p^ parameter) and accordingly 

have similar identities. 


The aerodynamic moment J (having the units of pound-feet and the 
dimensions L^MT - ^) and the aerodynamic pressure p (having the units 
of pounds per square foot and the dimensions L _ ^MP"^) are dependent upon 
the same physical quantities as the aerodynamic force R. Therefore, the 
functional relationship of the dimensionless parameters involving the 
aerodynamic moment J and that of the aerodynamic pressure p may be 
obtained by dimensional analyses similar to that for the aerodynamic 
force R, or, more conveniently, from equation (54) . Solving equation (54) 
for R yields : 



Moo 

W' 


9 



0oo4oo 




where the functional factor is the aerodynamic force coefficient. Since 
an aerodynamic moment is the product of a force and a lever arm 


J = ZR 


whence 



v«p l’ 

00 


, 


7 , 


> 


ejv, 





(55) 


Similarly, an aerodynamic pressure p is 
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so that 



Trfl’ 

00 


CpVrfi^l 


7 , 


®V^co 9oo4so 
3 ’ 5 J a ' 

v- pJ v. 3 Pc# i 




( 56 ) 


Similitude Equations 


With the model quantities denoted by the subscript m, and the air- 
craft quantities by the subscript a, the aerodynamic force, moment, and 
pressure similitude equations are obtained from equations (54), (55), 

-and (56) as: 



CpVcoPcpA _ f PpV-P-A 

\ ^ / m \ ^ /a 



( 57 ) 


Continued on next page 
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(57) 

Cone. 


Aerodynamic Force, Moment, and Pressure Relations Under the 
Conditions of Aerodynamic-Heating Similitude 

Equations (57) express the conditions of similitude that exist between 
an aircraft and a model with respect to aerodynamic forces, moments, and 
pressures. It is now desired to determine whether these similitude con- 
ditions are compatible with the similitude conditions of part IV that 
must be fulfilled in order to subject the model to heating by aerodynamic 
means in simulation of the aerodynamic heating experienced by the aircraft. 

The last two similitude equations of equations ( 57 ) are observed to 
involve the angular velocity & and the angular acceleration a. As dis- 
cussed in part XV, the aerodynamic method of heating the model is not 
capable of satisfying the heating similitude with respect to a, and a 
and accordingly provides heating simulation only when the portion of the 
heating attributable to a and co is negligible relative to the remainder 
of the heating. Consequently, a and ex are assumed to be zero, where- 
upon the last two similitude equations of equations ( 57 ) vanish. 

It is next observed that of equations ( 57 ) the fourth through the 
tenth are identical to the second through the eighth of the aerodynamic- 
heating similitude equations (40) . Therefore, with respect to these 
similitude equations, the imposing upon the model of heating by the aero- 
dynamic means of part IV simultaneously satisfies these similitude con- 
ditions with respect to the attendant aerodynamic forces, moments, and 
pressures imposed on the model. Accordingly, there remains only the 
first three of equations ( 57 ) "to be considered, which are 
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(58) 



(59) 



( 60 ) 


By the substitution of equations (5), (46), and (47), equations ( 58 ), 
( 59 ) > and ( 80 ) reduce to 


% = nRa 

( 61 ) 

11 

% 

** 

( 62 ) 

1 

% “ a Pa 

( 63 ) 


When the model is subjected to heating by the aerodynamic means of 
part XV, the forces, moments, and pressures experienced by the model are 
given, respectively, by equations ( 6 l), ( 62 ), and ( 63 ). These forces, 
moments, and pressures are to be compared with the forces, moments, and 
pressures found in part III to be required for simulation of internal 
stresses due to externally applied aerodynamic loade. 

For simulation of internal stresses due to externally applied aero- 
dynamic forces, the required force similitude is given by equation ( 32 ) 
of part III as 


F m = n2jr a 

From this equation and equation (6l), since R a = F a , it is found 
by substitution that 



( 64 ) 
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Thus, -when the model is heated by the aerodynamic means described in 
part XV, the aerodynamic forces imposed upon the model are l/n times 
as great as those required for simulation of internal stresses due to 
externally applied aerodynamic loads. 

For simulation of internal stresses due to externally applied aero- 
dynamic moments the required moment similitude is given by equation (54) 
of part III as 


N m = 

If this equation' is compared with equation ( 62 ), since J a = N a , it is 
found by substitution that: 


Thus, when the model is heated by the aerodynamic means discussed in 
part IV, the aerodynamic moments imposed upon the model are l/n times 
as great as those required for simulation of internal stresses due to 
externally applied aerodynamic moments. 

For simulation of internal stresses due to externally applied aero- 
dynamic pressures the required pressure similitude is given by equation (53) 
of part III as 


^m “ -^a 


Comparison of this equation with equation ( 65 ) shows that, since p a = F a , 


p m = 



(66) 


Therefore, when the model is heated by the aerodynamic means described 
in part IV, the aerodynamic pressures imposed upon the model are l/n 
times as great as those required for simulation of internal stresses due 
to exter nal ly applied aerodynamic pressures. 

Inspection of equations (64), ( 65 ), and ( 66 ) shows that, when n < 1, 
the forces, moments, and pressures imposed on the model will exceed those 
required for stress due to external aerodynamic loads, even when the cor- 
rections for inertial effects and gravity are disregarded. It is apparent, 
then, that the model subjected to heating in a wind tunnel will require 
that corrective loads be applied to it by some mechanical means. It 
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follows that the possib ili ty of testing the model for which n < 1 
in free flight, as discussed in part IV, is precluded by inability to 
correct the external loads to those required. 

Consider next the case where n = 1. This degenerate case is in 
essence the testing of the aircraft structure itself, and hence in the 
wind tunnel similitude to flight exists except for inertial loads and 
the effects of angular velocity and angular acceleration. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 15 , 195T- 
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